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I. INTRODUCTION 

Two-photon decay processes in hydrogen-like ions repre- 
sent an intriguing physical phenomenon and are the subject 
of intense research. The metastability of the 25 level, which 
is limited only by two-photon decay, makes it amenable to 
high-precision measurements. Interestingly, though, the two- 
photon decay has never been studied within the so-called Za- 
expansion beyond leading order, that is, beyond the order of 
a 2 (Za) e for the decay width in units of the electron rest mass 
(in this article, we use natural units, h = c = e = 1). 

The first study of the two-photon decay rate T of the 2S 
state was carried out by Goppert-Mayer in 1931 [1], and the 
well-known nonrelativistic result was derived, 

t- 1 = T Q = 8.229 352 Zh' 1 = 1.309 742 Z 6 Hz . (1) 

This result has been verified experimentally [2-4]. 

In the non-recoil limit, the leading correction terms mod- 
ifying this result are given by a relativistic correction of rel- 
ative order (Za) 2 and a radiative correction of relative order 
a (Za) 2 ln[(Za)~ 2 ]. We can write the following expansion, 
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with coefficients 72 and 73 to be determined. 

The next higher-order term not included in Eq. (2) is a 
nonlogarithmic radiative correction of order a (Za) 2 . Equa- 
tion (2) is complete up to order a 3 ln(a). 

The coefficient 73 is known for the 2S-1S transition [5, 6], 
but it remains unknown for any other two-photon transition in 
a hydrogenlike ionic system. The coefficient 72, which intu- 
itively could be assumed to represent an easy computational 
task, has not yet been calculated for any two-photon transi- 
tion, to the best of our knowledge. We address both 72 and 73 
in this paper. 

The relativistic correction of relative order (Za) 2 actually 
involves quite a large number of individual contributions: (i) 
multipole (quadrupole radiation) correction, (ii) relativistic 



corrections to the electron's transition current, and (iii) rela- 
tivistic corrections to the Hamiltonian and to the bound-state 
energies of initial and final states, due to zitterbewegung, rel- 
ativistic kinetic energy, and spin-orbit coupling. Each one of 
these contributions entails a computationally demanding sum 
over virtual states, and an integration over the photon energy. 
We here calculate the corrections one after the other and check 
gauge invariance all along the way. Finally, we obtain rigor- 
ous results for 72 and 73. 

Our approach is inspired by Nonrelativistic Quantum Elec- 
trodynamics (NRQED), albeit in a restricted way: in a two- 
photon decay, the photon energies are bound by the energy 
difference of the initial and final states, and therefore the prob- 
lem of separating the energy scales of the high-energy vertex 
terms does not arise. However, the interaction Hamiltonian 
still has to be expanded in the sense of NRQED, and we have 
the choice between two gauges which determine the form of 
the interaction Hamiltonian. Either the "length" (Yennie) or 
"velocity" (Coulomb) gauge can be chosen. The final result 
should not depend on the gauge. 

In the Appendix of Ref. [7], the gauge invariance of the 
two-photon decay rate was shown to hold within the fully rel- 
ativistic formalism, within the class of fully relativistic gauge 
transformations given by Eq. (A8) of Ref. [7]. The Power- 
Zienau gauge transformation [8] as given in Eqs. (18) and (19) 
of Ref. [9] has a nontrivial dependence on the coordinates and 
allows us to express the QED interaction Hamiltonian exclu- 
sively in terms of observable field strengths which in turn cor- 
respond to derivatives of the vector potential. This transfor- 
mation is most suitable for a nonrelativistic treatment, but due 
to the nontrivial dependence on the coordinates and due to 
problems related to the physical interpretation of non-gauge 
invariant quantities [10-12], a few subtleties arise. 

After considerable discussion on this point within the 
community [10-12], the conclusion has been reached that 
gauge transformations have to be considered very carefully 
in bound-state problems. E.g., for the radiative corrections to 
the two-photon decay rate [5], the results are invariant under 
a "hybrid" gauge transformation [11], where the interaction 
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Hamiltonian is gauge transformed, but the gauge transforma- 
tion of the wave function is neglected. I.e., although a gauge 
transformation normally entails a local, "pointwise" transfor- 
mation of the wave function, this whole transformation is 
flatly ignored, and the "usual" Schrodinger eigenstates [13] 
are used for initial and final states of the process under in- 
vestigation. We show here that the relativistic corrections to 
the two-photon decay rate are invariant under such a trans- 
formation (the gauge invariance of the leading logarithmic 
QED corrections was shown in Ref. [5]). In general, prop- 
erties of atomic states which can be formulated using adia- 
batic 5-matrix theory are invariant under this kind of hybrid 
gauge transformation, whereas in time-dependent problems, 
the choice of gauge has to be taken into account even more 
carefully [10-12]. In the latter case, the gauge transformation 
of the wave function cannot be ignored. 

When generalizing the results to higher excited initial and 
final states, one has to overcome a few subtle difficulties, be- 
cause one has to separate the 35-15 double-dipole (E1E1) 
two-photon decay from the cascade 35-2P-15. The 2P state 
appears both as a virtual state for the two-photon decay pro- 
cess as well as an intermediate state for the cascade process. 
In the two-photon decay rate, when regarded as differential 
with respect to the photon energy, the presence of the 2P state 
causes a (quadratic) singularity. Because we are interested in 
the total decay rate, we have to integrate over this singularity, 
which is quadratic and thus a priori not integrable. Remov- 
ing the 2P state from the sum over virtual intermediate states 
leads to gauge-dependent results [14—18]. In order to separate 
the cascade contribution from the two-photon correction for 
the two-photon decay, one has to use a special integration pre- 
scription detailed in Refs. [17, 19-21]; the prescription con- 
stitutes a generalization of the principal-value integration to 
quadratic singularities. Here, we extend the relativistic calcu- 
lations for two-photon decays to highly excited initial states 
using this formalism. 

We organize the paper as follows. In Sec. II, we explain 
the theoretical methods used in our approach. In Sec. Ill, we 
consider all the corrections separated by their physical origin 
for the 25-15 transition and show explicitly that each con- 
tribution is gauge invariant. In Sec. IV, we present numer- 
ical results for the 25-15 transition and also for transitions 
from higher excited states, and we discuss the separation of 
the cascade contribution from the coherent two-photon cor- 
rection to the decay rate. Results for the QED radiative cor- 
rections of logarithmic order are presented in Sec. V. Con- 
clusions are drawn in Sec. VI. As already mentioned, natural 
units h = e = c = 1 are used throughout this paper. 



II. THEORETICAL BACKGROUND 

The two-photon decay rate is given as the imaginary part of 
the two-loop self-energy correction [22] which can be derived 
using nonrelativistic Quantum Electrodynamics (NRQED) 
[23]. A detailed derivation of the nonrelativistic two-photon 
decay rate, valid for all transitions including those involv- 
ing highly excited states, is contained in previous works 



[5, 17, 19, 21], and there is no need to reproduce it here. 

We recall that in velocity (Coulomb) gauge, the interaction 
Hamiltonian for the interaction of the electron with the quan- 
tized radiation field is given as 



Hr = 
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p - A + A-p 



e 2 A 2 
2m ' 



(3) 



where p is the electron momentum, A is the vector potential, 
and m is the electron mass. This interaction leads to the fol- 
lowing expression for the nonrelativistic decay rate, 
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Here, "Re" denotes the real part, and the limit e — > is taken 
after all integrations have been performed. The summation 
convention is used throughout this article. The superscript £ 
denotes the velocity-gauge form of the expression. 

For length (Yennie) gauge, the (leading) interaction Hamil- 
tonian takes the simple form 



Hr = -eE-f. 



(5) 



If this Hamiltonian is used, we obtain for the nonrelativistic 
expression 
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Using the relation [10, 24] 
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the equivalence of these two expressions can be shown. Note 
that this is only valid if a complete spectrum is used for the 
representation of the propagator. 

For a fully relativistic calculation of the effect, we would 
have to use the Dirac Hamiltonian Hd in the propagators in- 
stead of the Schrodinger Hamiltonian H, and also the inter- 
action Hamiltonian and the wavefunction would have to be 
changed accordingly. However, as we want to work nonrel- 
ativistically, we transform the fully relativistic Dirac Hamil- 
tonian and its interaction Hamiltonian into effective nonrel- 
ativistic operators. This can be achieved by using a Foldy- 
Wouthuysen transformation [25], which identifies the nonrel- 
ativistic Hamiltonian as the leading term, and thus leads to a 
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systematic way of expressing the relativistic corrections. Fur- 
thermore, it allows us to express the relativistic corrections to 
the electron's transition current within the Za expansion. 

Alternatively, one can resort to the literature [26], where the 
corrections to the Schrodinger Hamiltonian have been tabu- 
lated. For the non-interacting part, this procedure leads to the 
well-known corrections to the Schrodinger Hamiltonian H, 
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The Darwin term proportional to the Dirac S originates from 
the zitterbewegung of the electron. The next term is the spin- 
orbit coupling, and the last is the correction due to the rela- 
tivistic kinetic energy. The relativistic corrections to the refer- 
ence state wavefunction and to its energy thus read as follows, 



E^E + SE = E+(<f> \SH\ $) 



E-H 



(9) 



5H\&) . (10) 



The transition current of the electron can be derived by acting 
with the Foldy-Wouthuysen transformation on a Dirac Hamil- 
tonian which is coupled to an electromagnetic vector poten- 
tial. The velocity-gauge result for the interaction Hamiltonian 
thus is (see Refs. [9, 27]) 
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We remember that the photon emission is characterized by the 
creation part of the electromagnetic vector potential operator, 
which carries a dependence of exp(— ifc • r). The transition 
current J can thus be written as 
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As we are considering a two-photon effect, contributions from 
seagull terms also have to be taken into account (here, two 
photons emerge from the same vertex). Terms proportional to 
A 2 are included in the seagull Hamiltonian which is given by 
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Expanding in powers of (Za) and extracting the photon cre- 
ation part, we obtain the seagull correction in relative order 

(Za) 2 , 



2m y 2m 3 Am 3 



(14) 



written in such a way that it multiplies the (creation part of 
the) photon fields A 1 AK 

The interaction Hamiltonian in length gauge, including rel- 
ativistic and multipole corrections, can be obtained by em- 
ploying two consecutive Power-Zienau transformations [8] 
after the Foldy-Wouthuysen transformation. This has been 
shown in Ref. [28]. The interaction Hamiltonian in length 
gauge thus reads 
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Here, the subscript separated by commas denotes the spatial 
derivatives with respect to the indicated Cartesian coordinates, 
evaluated at the origin [28], which is defined to be the loca- 
tion of the ionic nucleus. This corresponds to a length-gauge 
transition current 
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where {A, B} — A B + B A is the anticommutator, and we 
examine the emission of a photon with four- vector (w,k). We 
are now in the position to discuss how the corrections to the 
decay rate can be determined from the transition currents in 
the two different gauges. We start with the velocity gauge. 

A. Velocity Gauge 

The nonrelativistic two-photon decay rate in velocity gauge 
[see Eq. (4)] can be written as 
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where the superscript £ denotes the velocity-gauge expression. 
Here, due to energy conservation, ui 2 = E^, i — E<p f — u)\, and 
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For the gauge invariance of this nonrelativistic expression, see 
Eq. (7). We only remark that the statement of gauge invariance 
can be brought into the compact from 



£ = -wi uj 2 C • 



(19) 
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where ( is defined in Eq. (26) below. Here and in the follow- 
ing, we suppress the superscripts ij of the £ and ( tensors in 
order to ensure the compactness of the notation, and we imply 
that £ 2 = £ y £ y (the indices i and j are summed over), and 
that = £ y S^ tJ . We define 5£ to denote the sum of the 
corrections due to all the previously discussed perturbations 
(Hamiltonian, energy, and current) and express the first-order 
relativistic correction 5T to the decay rate as (see Ref. [5]) 

4a 2 f 4c? f 
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The correction Suj max = SE<s> i — 5E<j, f is necessary to ensure 
that the perturbed energy conservation condition is fulfilled: 

uji + uj 2 = E$ i - EiSfj + 5u) max , (21a) 

= <*i|«ff|*i> - (*f\5H\* f ) . (21b) 

so that the frequencies of the two quanta add up to the per- 
turbed transition frequency. However, due to the presence of 
the seagull terms, further corrections have to be taken into ac- 
count. 

After some algebra, we see that <5£ can be expressed as the 
sum of fifteen terms that account for all the relativistic and 
multipole perturbations, 



5£ = £>&. 



(22) 



fe=i 



The perturbations of the energies of the initial and final states 
lead to the following terms, 
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The perturbations to the initial and final-state wavefunctions 
lead to the following four effects, 
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terms (observe the different denominators), 
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The correction to the electron's transition current can affect 
both the initial and the final states, and this gives rise to a total 
of four terms, 
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The seagull Hamiltonian acting on the unperturbed wavefunc- 
tions leads to 



5Z 13 = -($ f \5S ij \$t) . 



(23m) 



The minus sign originates because we have written all matrix 
elements (second-order perturbations) in the "1/(H — E)" 
form, which corresponds to a negative second-order energy 
perturbation. In order to be consistent, we have to use the 
negative higher-order seagull Hamiltonian, which is applied 
in first-order perturbation theory. Finally, we have the seagull 
terms which were already present in Ref. [5] which account 
for the emission of two photons from the perturbed initial state 
or to the perturbed final state. They are given as 
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where we invoke second-order perturbation theory with the 
leading seagull term e 2 A 2 /(2m). Using a complete basis- 
set of hydrogen eigenfunctions and their orthonormality rela- 
tions, we can show that 



(24) 



The reason is that both (5£ 14 and <5£i 5 are proportional to the 
non-diagonal matrix element (<&f\5H\$i), but with opposite 
pref actors. 
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B. Length Gauge 

The nonrelativistic, length gauge expression in Eq. (6) can 
be written as 

E$>^ — E$> j. 

rC = ^7 / ^i^2 3 C 2 , (25) 
o 

where the superscript ( denotes the length-gauge expression. 
Here, uj 2 is defined as in Eq. (17), and 
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Following the same procedure as for the velocity gauge ex- 
pression, we can write the first-order correction to the two- 
photon decay rate in length gauge, 
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(27) 

where again 5( denotes the sum of all the correction terms in- 
curred by the relativistic perturbations of the Hamiltonian, and 
of the energies of the initial and final states, and of the length- 
gauge current. Indeed, in the length gauge, the correction S( 
contains only twelve as opposed to fifteen terms, 
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The energies of the initial and final states are perturbed and 
this gives rise to the first two correction terms, 
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In complete analogy to Eqs. (23c) — (23f), the perturbations to 
the initial and final state wavefunctions are accounted for by 
the following four terms, 
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Furthermore, the corrections from the perturbed Hamiltonian 
give rise to two terms, 
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The length-gauge correction to the current SI gives rise to four 
more terms, 
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The seagull term is not present in the length gauge. In the next 
section, we analyze these corrections in the light of gauge in- 
variance. We separate the corrections by their physical origin, 
and show more than the gauge invariance of the final result: 
namely, we are able to demonstrate that each physically dis- 
tinguished correction is gauge invariant in itself. 



III. GENERAL PROOF OF GAUGE INVARIANCE 
A. Orientation 

First of all, let us remember that in all bound-state calcu- 
lations, we actually use a hybrid gauge transformation [11] 
where we ignore the gauge transformation of the wave func- 
tion. The non-interacting relativistic Hamiltonian, given in 
Eq. (8), by definition is gauge invariant. Thus, we only gauge 
transform the electron's transition current and the photon field 
operator, or alternatively, we let the interaction Hamiltonian 
undergo a gauge transformation. We show here that the full 
gauge invariance is obtained by carefully considering the in- 
terplay of the relativistic corrections to the wavefunction, to 
the Hamiltonian and to the energies of the bound states (the 
initial and the final states). 

The whole problem becomes simpler when it is divided 
into three distinct parts, the first of which is a generalized 
correction due to the relativistic Hamiltonian, the second of 
which is a quadrupole correction, and the third is a remaining 
correction (a further correction to the current), which can be 
shown to vanish after the use of commutator relations. Gauge 
invariance can be shown for each of these corrections sepa- 
rately, provided some parts of the velocity gauge correction 
to the electron's transition current (12) are identified as be- 
ing generated by the relativistic Hamiltonian (8), and treated 
together with the correction to the Hamiltonian. Here, the 
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velocity-gauge expression appears to be more complicated. 
The quadrupole correction, by contrast, looks a little more in- 
volved in the length gauge. Gauge invariance with respect to 
the velocity gauge can be shown provided we include a part 
of the seagull term (14) into the velocity-gauge expression for 
the quadrupole term. It is then relatively easy to show that all 
remaining terms vanish separately. 

In the following, we discuss the general approach to the 
proof of gauge invariance in some detail. Further aspects are 
elucidated in Appendices A and B. 



B. Correction to the Hamiltonian 

Let us discuss first the general paradigm and start with the 
corrections induced by the relativistic Hamiltonian (8). The 
gauge invariance for the leading-order term (the nonrelativis- 
tic result) can be traced to the formula (19), 



-OJl LU2 C , 



(30) 



where £ represents the velocity-gauge form and £ represents 
the length-gauge form. 

Both £ and £ actually carry superscripts ij , which we sup- 
press here to leave the notation compact, as already discussed. 
Let us now suppose that the total velocity-gauge correction 
due to the relativistic Hamiltonian can be expressed as 8£h, 
and the corresponding length-gauge expression is 5£h ■ The 
precise definition of S^h and S(h will be discussed later. We 
are able to show the following gauge invariance relation, 



(31) 



based on which we can prove the gauge invariance of the en- 
tire correction 5Th due to the relativistic Hamiltonian, 
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Here, again, the superscript £ denotes the velocity gauge, 
whereas £ denotes the length gauge. We are indeed able to 
show such a relation for all three terms given in (8), but only 
if we include in the definition of S^h specific corrections to 
the electron's transition current. Our gauge-invariance rela- 
tion can be illustrated as follows. The correction <5£# con- 
tains the wavefunction correction in the velocity gauge, the 



Hamiltonian correction in velocity gauge, the energy correc- 
tion in velocity gauge, and the seagull term in velocity gauge, 
as well as the current correction due to the current operator 
5J l H = — i[r\SH}. By contrast, 5(h equals the sum of the 
wavefunction correction in length gauge, the Hamiltonian cor- 
rection in length gauge, and the energy correction in length 
gauge. Note that the term — Scj max uji ( in Eq. (31) is related 
to the modified energy conservation condition, and that Suj max 
here is the correction to the transition frequency due to the rel- 
ativistic Hamiltonian given in Eq. (21b). Using this result, we 
are able to show that the total correction to the decay rate due 
to all three terms given in (8) is gauge invariant. 
The current that we add in the velocity gauge is 
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The seagull term that we add in velocity gauge is due to a 
double commutator 
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This term is part of the seagull Hamiltonian (14). We are now 
in the position to give the precise definition of 5£h and S(h, 
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For further details, see Appendix A. 



C. Quadrupole (multipole) correction 

The quadrupole correction is not associated with any cor- 
rection to the bound-state energy or to the Schrodinger Hamil- 
tonian. It can be treated separately and identified with a cor- 
rection 6 Jq to the current in velocity gauge, and with a cor- 
rection SIq in length gauge. The velocity-gauge current is 



1 P i n -A2 
- — (k ■ r) . 

2m l ; 



(37) 



We can ignore the first term because it vanishes after angular 
algebra, for the first-order correction to the two-photon decay. 
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This is unlike the (Za) 2 correction to the Lamb shift, where 
this term contributes as a simultaneous perturbation to both 
currents, because one and the same photon is being emitted. 
Here, two photons are being emitted, and angular averaging 
occurs for both of them separately. 

The quadrupole current in the length gauge is 



51} 



-k- 



1 



1 



2muj 



(l x ky 



6niLu 



(Lxk) l (k-r) + (k-r)(Lxk) 1 



-- 6 (k- r r 



6niLU 



(L x fc)*(fc -r) + (k- f)(L x kf 



(38) 



where in the last step we have ignored the terms that vanish 
after angular integration. We find that the quadrupole term 
is gauge invariant provided we include, in the velocity-gauge 
expression, the seagull contribution from the term 



8S% 



- — (k-f) 2 5 ij . 
2m v ' 



(39) 



Now, the sum of SSq and SS 1 ^ is the full higher-order seagull 
term SS tJ given in Eq. (14). 

We denote the correction to the quadrupole matrix element 
in the velocity gauge by 5£q (it includes the seagull correction 
due to SSq) and use 5(q for the corresponding correction to 
the matrix element in the length gauge. We are able to show 
that 



Eis, —En> 



Aa z f 
2— / dwiWiw 2 £(f£Q 

o 

Eq ■ — Eq j. 

4a 2 ' f 

2 ~^T / ^i^iw 2 (-wiW2C)[- w i^2'5Cq] (40) 



Eq^ — Eq j. 



proving the gauge invariance of the quadrupole correction. 
The precise definition of <5£q and S(q reads as follows, 



12 



13 



SJ=SJ Q 



(41) 



8S=8Sq 



and 



i=9 



(42) 



6I=6I Q 



Further details are provided in Appendix B. 



D. Remaining corrections 

We have by now treated the correction due to the en- 
tire Hamiltonian (8), the entire seagull term (14) and the 
quadrupole interaction. The remaining terms are current cor- 
rections. In the velocity gauge, these read 



sji = sr 



5Jh 



SJ, 



Q 



= -2^(* X *) , ( 1 " i *- i ')- 



1 Za 

Am 2 r 3 



(fx a) 1 . 
(43) 



Using commutator relations, it is possible to show that 
1 



P 



H-E^.+w 



$1 



(a x k) 3 



(a x ky 



1 



4\; 



0. 



(44) 



This relation is valid for both k = k\ t 2 if w is changed accord- 
ing to Eqs. (23i) to (231), and for an arbitrary initial and final 
state. Thus, the contribution of the first term on the right-hand 
side of (43) vanishes. Furthermore, we can replace 

-i(" xfc T(H^-4^ x ^' 



1 Za i 

(r x a) 



Am 2 



1U) 

Am 2 



(<? x v) 1 , 



(45) 



when contracted with the photon propagator. This relation is 
known from Lamb shift calculations (see Ref. [29]). There- 
fore, the entire contribution from the remaining corrections to 
the current vanishes in the velocity gauge. 

In the length gauge, the remaining corrections to the current 
are given as 



Sli = (a x jfeY ( 1 - ijfe • f) + — (a x r) 1 



(46) 



The first term vanishes in view of Eq. (44). The remaining 
terms also do not contribute to the corrections to the decay 
rate. This follows from the relation 



(47) 



for the last two terms of Eq. (46) when contracted with the 
photon propagator. The precise definition of <5£r and 5(r 
reads as follows, 



12 



(48) 



8J=SJ R 



and 



12 



(49) 



SI=SI R 
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IV. NUMERICAL CALCULATIONS 

A. 2S-1S Decay 

The phenomenologically most important two-photon decay 
process is the 2S-1S decay. Our gauge-invariant result for the 
correction to the decay rate due to the relativistic Hamiltonian, 
as discussed in Sec. IIIB, reads 



ST H = To [-0.5082 {Zaf 



(50) 



For the quadrupole correction, the gauge-invariant result is 
(see Sec. IIIC) 



ST Q = T [-0.1555 {Zaf 



(51) 



The remaining current corrections vanish, as discussed in 
Sec. HID, 



STr = . 



(52) 



The total result for the relativistic correction to the two-photon 
decay rate thus reads 

ST = ST H + ST Q + ST R = T [-0.6636 [Zaf] . (53) 

It is instructive to break down the corrections to the Hamilto- 
nian further. Namely, according to Eq. (8), we have the zitter- 
bewegung (zb) term, 



nZa 3 



2m 



the kinetic energy (ke) term, 



8m 3 ' 

and the spin-orbit (LS) coupling 

SH LS = Za 

Am 2 r 3 

The corresponding results read, for the 25-15 decay, 

<5r zb = r [-0.7577 {Zaf] , 
<5r ko - r [0.2495 {Zaf} , 

<5r LS = o. 



(54) 



(55) 



(56) 



(57a) 
(57b) 
(57c) 



This concludes our discussion of the two-photon decay of the 
2S state, and we can now proceed to calculate decays from 
higher excited states. 



B. Higher Excited States 

In principle, one might assume that in order to calculate the 
relativistic correction to the two-photon decay from higher ex- 
cited states, only the initial and final state wavefunctions have 



TABLE I: Results for the 72 coefficient as defined in Eq. (2). This 
coefficients gives the relativistic corrections to the two-photon decay 
rate. 







|*/>=|lSl/2> 


l*/> = |25 1/2 > 


I*i> = 


I2S1/2) 


-0.6636 






I3S1/2} 


-2.6637 


-1.7038 


1*0 = 


|4Si/ 2 > 


-4.5192 


-7.8530 


1*0 = 


|3£> 3 /2> 


-2.2978 


7.8533 


1*0 = 


1 3^5/2 > 


-1.0981 


-22.2671 



to be changed accordingly. However, historically the gener- 
alization to higher excited states has proven to be problem- 
atic. For higher excited states, the two-photon transition can 
take place not only through virtual intermediate states with an 
equal or higher energy than the initial state, but also through 
cascades via intermediates states with a lower energy. For the 
35 initial state, a decay via the cascade 3S-2P-1S is possi- 
ble. The allowed cascade transitions cause singularities in the 
propagators. As we are interested in the total decay rate, we 
integrate over the propagators and thereby also over the sin- 
gularities. These singularities are quadratic and thus a priori 
not integrable. 

Finally, after some discussion [14—19, 30], the conclusion 
has been reached that the two-photon correction to the decay 
width of the initial state can be obtained using an integration 
prescription where the double poles are treated in a manner 
inspired by quantum electrodynamics, where the photon en- 
ergy integration contour extends infinitesimally into the com- 
plex plane [21, 31]. Note that the two-photon correction thus 
obtained is a further correction that has to be added to the 
one-photon decay width that is otherwise responsible for the 
cascade transition. Using this procedure, we were able to de- 
termine the relativistic and multipole corrections to the non- 
relativistic decay rate for many higher excited states which 
fulfill the same gauge relations as for the 25-15 transition. 
Final results are given in Table I. 



V. LEADING LOGARITHMIC QED CORRECTIONS 

The zitterbewegung term in the relativistic Hamiltonian, ac- 
cording to Eq. (54), is given as 5H zh = nZa 5 3 {r)/{2m). 
The effective potential that gives the leading QED radiative 
corrections is 



<5-ffrad = 



4a 



{Za) \n[{Za)- 



8 3 {r) 



(58) 



This relation implies that the 73 coefficient can be obtained as 
8 72,zb/3 where 72, z b is the contribution to 72 caused exclu- 
sively by the zitterbewegung term. As this contains no spin 
dependence, the 73 coefficient is spin independent. For the 
25-15 transition, e.g., we have according to Eq. (57a), the 
relation 73 = §(-0.7577) = -2.0205. Results for other 
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transitions are given in Table II. The 73 coefficient becomes 
numerically rather large for 3D-2S transitions. Note that the 
correction is the same for decay from 3D 3 / 2 and 3D5/2 be- 
cause the potential (58) does not involve any spin-dependent 
terms. 



TABLE II: Results for 73 as defined in Eq. (2). 







1*/) = \is 1/2 


) 1%) = \2S 1/a ) 


I*i> = 


1 251/2 ) 


-2.0203 




l*i> = 


j 35*1/2 ) 


9.6521 


16.0424 


l*i> = 


|45 1/2 > 


20.7364 


61.7499 


l*i> = 


1 3^3/2) 


-5.4681 


144.3639 


l*i> = 


|3 J D B /2> 


-5.4681 


144.3639 



VI. CONCLUSIONS 

The precise treatment of the two-photon decay width in 
ionic hydrogenlike bound systems with low nuclear charge 
numbers demands an evaluation of the relativistic and mul- 
tipole correction of relative order (Za) 2 , which is the lead- 
ing correction to the classic result [1]. The leading logarith- 
mic QED correction of relative order a(Za) 2 \n[(Za)~ 2 ] 
also needs to be determined. These corrections can be 
parametrized according to Eq. (2) in terms of two coefficients 
72 and 73 which are given in Tables I and II. 

Of particular interest is the result 

72(25- 15) = -0.6636 (59) 

for the 25-15 decay. This result [see Eq. (53)] is the sum 
of a correction due to the relativistic Hamiltonian [Eq. (50)] 
and a correction due to the quadrupole term [Eq. (51)]. We 
also generalize our approach to the two-photon decay from 



higher excited states (Tables I and II). As usual in quantum 
electrodynamic calculations, the magnitude of the correction 
terms grows with the principal quantum number. The decay 
from D states is also treated, and it is worthwhile noting that 
the spin-independent logarithmic correction terms of relative 
order Za 2 \n(Za) turn out to be large in magnitude (see Ta- 
ble II). Finally, as shown in Appendix C below, a comparison 
of our results to those of a nonperturbative (in Za) calcula- 
tion for the 35-15 decay (Ref. [20]) reveals that the term of 
relative order (Za) 2 can account for the bulk of the relativis- 
tic correction up to some rather high nuclear charge numbers 
CZ<40). 

With our NRQED-inspired approach, we can uniquely 
identify the physical origin of the (Za) 2 -correction terms to 
the two-photon decay width, as discussed in Sees. 111 B, IIIC 
and HID, and give their values separately. It is sometimes 
worthwhile to use the effective nonrelativistic treatment of 
NRQED, because it may yield information which could not be 
obtained by a fully relativistic treatment, regarding the break- 
down of the corrections. Furthermore, the calculation of the 
full spectrum of the propagator can be greatly simplified us- 
ing lattice methods [32], increasing the speed as well as the 
numerical stability of the evaluation, which is especially im- 
portant in the domain of low nuclear charge numbers. 

Another aspect is that the proof of the gauge invariance, as 
carried out in full detail in Appendices A and B, turns out to be 
a surprisingly lengthy calculation. We stress once more that 
the gauge invariance is shown to hold even if we ignore the 
gauge transformation of the wave function, in the sense of the 
"hybrid" gauge transformation developed in Refs. [11, 12]. 
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I 

APPENDIX A: GAUGE INVARIANCE OF THE HAMILTONIAN CORRECTION 



We give further details regarding the gauge invariance of the seagull term. Useful general relations are p 4 = im [H — E + ui , r l ] 
and ui2 = E<s, i — E& f — u\. The term 6£i can be transformed to 



^1 



J \V 

m \H — E$ i + UJ\ J m 
1 



$ < )<$ j |<yff|$ i > = -a;iW2(*/ 



+ (u>2 - Wi) < $/ 



H - E 9i + uji 
An analogous relation also holds for 5£ 2 , 



H - E^ + wi 
) (QilSH^i) + ($ f |rV| $i> ($i|Aff|*i) . 



1 



m \H — E&, — u>i J m 
1 



$i } ($f\6H\$f) = -ll>ilu 2 ( $/ 



+ (wi -w 2 ) ( 



H — E<s> f — uji 



H — E<£ f — ui\ 
) ($ f \6H\if> f ) + ($/ |rV|*i> (& f \6H\$ f ) . 



(Al) 

*^{* f \5H\* f ) 

(A2) 
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These relations are equal to those found in Ref. [5] for a radiative correction potential. The relations for the correction to the 
wavefunctions are altered because we are considering a different Hamiltonian. Thus, ££3 gives 



m H 



— E<s> i + u>i m \ E,j> i — H 



SH 



r Ex 



- W2 ( ®f 



=T 3 



- |rV| ($j \SH\ + <$/ |rV'5H| $ 4 ) + ^2 
For 5^4 this yields 



r^SH 



(A3) 



^4= (*/ 



m 



H - E 9f - u>i m \E^ t - H 



SH 



' I r j ( ) SH 



— Etp f — u>i 



SH 



$4 



=T 4 



For the correction <5£5 to the final-state wavefunction we get 



^(,22 



(A4) 



2?$ f — 22 / m 22 — E^ i + ui m 
1 



$^<$, |<522| $/> + /$/ 



57? 



5221 



i r t i , 

^E^j—HJ H — E<$ i -\-u)\ 



22 - E^ i + uji 
- ($/ \SH\ *,) ($/ |rV| + ($/ |<JJfrV| - wi 



=T 5 
1 



22 - Eq,, + U>! 



(A5) 



and for <5£ 6 



SH 



— t — y_ 1 

E$ f — H J m H — E^ f — uii 



EL 
m 



SH 



+ u 2 (<bf\5H\* f )U f 



I 

H — — ui 



— — l r r 

Efyj—HJ H — E^j—LUi 



4>, 



=T 6 



However, 

transformed to give 



{<i> f \SH\^ f ){^ f \r l rn^) + {^ f \SHr l ^\^)-LU 2 (<i> f , , . . 

\ " — — w i / 

the corrections to the wavefunctions lead to some remainder terms which have to be analyzed separately, 
ed to give 



(A6) 
'. They can be 



T 3 +T 4 = =-(*/ 
m \ 



T 5 +T 6 = = — ( 



5i7 



- ($/ |rV| $,) ($i $,) - ($/ |rV J '<5i7| $,) , 



(A7) 



SH 



E<f>, — H 



$ 4 ) S ij + ($/ |<5iJ| $/) <$/ |rV J '| $j) - ($/ |jfTrV| *i) . (A8) 
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We observe the seagull terms <5£i4 and <5£i5 emerge and cancel, explicitly. The other terms on the right-hand side will be treated 
separately, later. The term ££7 arising from the correction of the Hamiltonian can be brought into length gauge form in the 
following way: 



p' 



-5H- 



V 3 



- lo 2 (®f 

Finally, for 8£s we nave 

p° 1 



m H - E^ i + lo\ H - E$ s Iwim 
' SHr^ 



<P t ) = lo y lo 2 ( $/ 



-6H- 



1 



H - + uji 



5£ 8 = - {®f 



-6H- 



r l 8H 



H-E^+LOi 



m H — E^ f — uj\ H — E<j> f — lo\ m 



Our intermediate result thus reads as follows, 



) +uj 2 { $/ 



r l SH 1 



H-E^+w-l H-E&1+GJ! 
$^ -{Qf^SHri^i) . 

1 1 

5H- 



(A9) 



H — E<$ f — o^i 



H — E^j, —u>i H — E^j—uj\ 



(A10) 



8 8 
i=l i=l 



5> 



C + L0 2 ($/ 



\6H] 



iJ — E^ — LOi 



4>, 



+ w 2 ( $/ 



-ff — E$j: — LOl 



$ i )-($ / |[[r i ,^],r^]|$ i ) . 

(All) 



where <5w max is defined in Eq. (21b). Fortunately, we can rewrite the terms with the [r 3 ,5H] commutators further, 



+ w 2 ( $/ 



H-E^+wi 
1 



1 p> 



if 



H - E$ ; + uji m 
f 1 



6J 



+ L0 2 { $/ 



1 



SJh 



H — E<&f — u>i 
1 



1 p> 



[r j ,8H] 



= -E*< 



i=9 



m H - E $i + Ui 

+ 2(<S> f \{{r\5H],ri}\<t> l ) . 



H — E<s> s — loi m 
p i 1 



(A12) 



m H — E^ f — LOi 



5,1 



H 



$ i ) + 2($ f \[[r i ,6H],ri]\$ i ) . 



6J=8.J H 



The current J# = — i[r J , SH] is defined in Eq. (33). Combining (Al 1) and (A12), we obtain the relation 



i=9 



E s & + E 5 & = ^^1^2 E ^* _ &mmWi c - 



6J=5J H 



i=l 



(A13) 



With the definitions (35) and (36), this leads directly to our gauge invariance relation (31). 



APPENDIX B: GAUGE INVARIANCE OF THE QUADRUPOLE CORRECTION 

For the proof of gauge invariance of the quadrupole correction it is more convenient to start from the length gauge expression. 
As the quadrupole term is a correction to the transition current, only the terms <5(9...i2 are relevant. The length-gauge transition 
current 61 is [see Eq. (38)] 



SPr 



/ bniLo 



(L x kfi-ik ■ r) + (-ijfe • r)(L x k) 1 



(Bl) 
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It is helpful to rewrite the second part of the transition current as 

(L x k)\-\k ■ r) + (-ifc • r)(L x kf = (fc ■ f)p\-\k ■ r) - r\k ■ p)(-ik ■ r) + (-ifc • r)(k ■ r)p l - (-ifc • f)r\k ■ p) . 
Using this and the general relations from Appendix A we can transform the first term <5( 9 to give 



(B2) 



-luiuj 2 5( 9 = Lo 2 - { $/ 



= $ 



p' 



1 



to H — + uji 



wi(jfei • r)V+ 3 \(L x fci) J '(fci • f) + (fci • r)(L x 
1 



TO 



m 



6 



r l (fci -r)V 



4\ 



— (ki-r) 2 ^ 

TO 



1 £> m 



:- L r l — r m — ^ri—r" 1 + rV m r'— — r"V m r J — 

TO 



TO 



TO 



TO 



(B3) 



For 5Cio we obtain in an analogous manner 



-wiw 2 <5Cio = ( <!>,• 



1 



cu 2 (k 2 ■ r)V+ - f(L x fc 2 ) J (fc 2 • rO + (fc 2 • x fc 2 ) j 



= (*/ 



m H — E<s> } — lui 



1 /- 



— (k 2 -r) 2 ^ 

TO 



gfe 2 ^ 2 



rV — r m — rV — r m + rV"V- ■ - /•'/•"'/"' — 



TO 



TO 



TO 



TO 



(B4) 



For the correction <5£ii with the current acting on the left side this yields 



(w 2 (fc 2 -f)V + /,■_,) '(7,2 r) - (A'2 ■ /" Mi. > /,- 2 )' 



1 /- 
2 



(£2 • r) 

r l (k 2 -r) 



1 



2^ 



i? — _E<j, i + wi to 
i 



1 



$») + Q^2[^>f r i (k 2 -r) 2 r j 



H-Eqi+U! 



: h "2 h "2 



jP_ r m r j _ r iP_ r m r j _|_ r ?n r lP_ r j _ r ?n r iP_ r j 



TO 



TO 



TO 



TO 



(B5) 



and finally for 5^12. 



-wiw 2 (5Ci 2 = w 2 g ( 



wi(Jfei • r)V+ - (Lx fci) J (fci • f) + (fci ■ f)(L x fci)' 



H — E^f 



r'(fci-r) 



2^ 

TO, 



i? — — u)\ to 
i 



r*(fci • r)V 



r l E- r m r 3 — r i ?— r m r i _|_ r m r lP_ r j _ r .m r »P_ r J 

mm mm 



(B6) 



Combining these results, we get 

12 12 

i=9 5I=8Iq , = 9 



SJ=SJq 



r V ( — r m — rV J — r m + rV"V — 

TO TO TO 



TO TO 



- ( ] ( <*>/ 



rV z — r m — rV J — r m + rV m r z — — rV" 1 ^- — r l r m r 3 

TO 



>P_ r m + r i r m r l ^_ 

TO TO 



TO TO 



+ ( 6 W ) ( */ 



TO 



r l E- r m r 3 — r i ?— r m r j + r m r l E- r i — r m r iP_ r 3 _ j-i^j-mP 3 
TO 



TO 



TO. 

z 



TO TO TO TO TO 



TO. 



(B7) 
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In order to simplify the resulting expression, we now commute the momentum operators in the remainder terms to the right, 



12 

i=9 



12 



SI=SIq j = g 



SJ=SJq 



(B8) 



The last two terms can be identified as the negative of the quadrupole contribution to the higher-order seagull term as given in 
Eq. (39), summed over the two photon momenta fci and fc 2 . Finally, this leads to the equality 



12 

—Ll>iLU 2 

i=9 



SI=5Iq i—g 



1:5 



SJ=SJq 



(B9) 



SS=SSq 



which verifies the gauge invariance relation given in Eq. (40). 



APPENDIX C: COMPARISON OF ANALYTIC AND 
NUMERICAL RESULTS 

We would like to compare our results for the analytic co- 
efficients listed in Tables I and II to numerical data obtained 
for 25-15 (see Ref. [7]) and 35-15 (see Ref. [20]). The au- 
thors of Ref. [7] obtained a fit to a convenient functional form 
in Za, leading to an approximate formula valid across the 
whole range of nuclear charge numbers Z [see Ref. [7] and 
alsoEq. (4.16) of Ref. [33]], 



To 



1 + 3.9448 (Za) 2 - 2.040 (Za) 4 
1 + 4.6019(Za) 2 



(CI) 



Upon re-expansion in Za, one may thus hope to obtain an 
estimate for the correction of relative order (Za) 2 . Indeed, 
the estimate thus obtained, 72 « —0.6571, is in fair agreement 
with the precise result (59), which reads 72 = —0.6636. 
For the 35-15 decay, we compare to a fully relativistic cal- 



culation carried out in Ref. [20], where the relativistic effects 
have been calculated for different values of Z. When using 
our results for 72, one can determine the corrected decay rate 
for different values of Z. For Z = 40, our analytic results aug- 
mented by the relativistic correction of relative order (Za) 2 , 
lead to a result of T « 1.61 (Z = 40) 6 rad/s, to be compared 
with the result T = 1.60 (Z = 40) 6 rad/s from Ref. [20] for 
the E1E1 two-photon decay rate. 

In general, there is quite a subtle interplay of the fully 
relativistic calculations with the Dirac-Coulomb propagator, 
which have meanwhile been done for a number of QED and 
other problems, and the Za-expansion approach: numerically 
more accurate results can be obtained with the former, and 
these are relevant especially for highly charged ions, but the 
physical origin of the relativistic corrections is much more 
transparent within the Za-expansion. Furthermore, the an- 
alytic calculations allow for a systematic expansion in powers 
of a and Za, as demonstrated in Eq. (2). 
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